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Abstract 

QCD plasma instabilities appear to play an important role in the equilibration of quark-gluon 
plasmas in heavy-ion collisions in the theoretical limit of weak coupling (i.e. asymptotically high 
energy). It is important to understand what non-linear physics eventually stops the exponential 
growth of unstable modes. It is already known that the initial growth of plasma instabilities in 
QCD closely parallels that in QED. However, once the unstable modes of the gauge-fields grow 
large enough for non-Abelian interactions between them to become important, one might guess that 
the dynamics of QCD plasma instabilities and QED plasma instabilities become very different. In 
this paper, we give suggestive arguments that non-Abelian self-interactions between the unstable 
modes are ineffective at stopping instability growth, and that the growing non-Abelian gauge fields 
become approximately Abelian after a certain stage in their growth. This in turn suggests that 
understanding the development of QCD plasma instabilities in the non-linear regime may have 
close parallels to similar processes in traditional plasma physics. We conjecture that the physics of 
collisionless plasma instabilities in SU(2) and SU(3) gauge theory becomes equivalent, respectively, 
to (i) traditional plasma physics, which is U(l) gauge theory, and (ii) plasma physics of U(l)xU(l) 
gauge theory. 



I. INTRODUCTION 



At high enough energy, heavy ion colhsions are presumed to create a quark-gluon plasma 
in local equilibrium. The physics associated with the creation and equilibration of this 
plasma is rather complicated, involving many different physical processes playing important 
roles at different stages of the collision. It is theoretically useful to try to understand the 
collision in the formal limit of arbitrarily high energies, for which the running coupling 
constant as is arbitrarily small due to asymptotic freedom. In this limit, it is believed that 
very early times are described by the saturation scenario [1-6], in which the non-equilibrium 
plasma starts out at early times tq ~ l/Qs as a gas of gluons with (i) non-perturbatively 
large phase-space density / ~ l/os and (ii) momenta of order a scale Qg known as the 
saturation scale. In the weak coupling limit as(Qs) "C 1, Baier, Mueller, Schiff and Son 
[7] attempted to find the parametric dependence of the equilibration time Tgq on CKg and 
found Teq ~ a~^^^^Q~^ . Their analysis is known as the bottom- up thermalization scenario, 
and the result arises from a complicated interplay of (i) the one-dimensional expansion of 
the plasma between the large, Lorentz-contracted, nuclear pancakes after the collision and 
(ii) a variety of individual coUisional processes which relax the plasma toward equilibrium. 
However, collective processes are often more important in plasmas than individual collisions. 
In particular, Mrowczyhski and others [8-18] have long suggested that plasma instabilities 
might play an important role in the equilibration of quark-gluon plasmas. Arnold, Lenaghan, 
and Moore [19] showed that this is indeed the case for the bottom-up analysis. That means 
that the bottom-up scenario of Baier et al. needs to be completely reanalyzed. The theory 
of quark-gluon plasma equilibration in heavy-ion collisions is currently in the embarrassing 
state of not even knowing how the equilibration time depends on in the weak coupling 
limit — that is, not even the power n is known in the parametric relation 

Teq ~ a-^Q:\ (1.1) 

For ultra- relativistic, homogeneous, parity-invariant, coUisionless plasmas, a long- 
wavelength magnetic instability known as the Weibel or filamentary instability generically 
arises whenever the momentum-distribTition of charged particles in the plasma is anisotropic. 
A proof (and a more precise statement) may be found in Ref. [19]. Reviews in the quark- 
gluon plasma literature of the physical origin of the instability may be found in Refs. [11, 19]. 
The actual quark-gluon plasma in heavy-ion collisions is not coUisionless, nor precisely ho- 
mogeneous, but Ref. [19] showed that the instabilities that plague the original bottom-up 
scenario are associated with small enough distance and time scales that collisions and in- 
homogeneity can be ignored.^ The analysis of such instabilities in the quark-gluon plasma 
literature has generally been restricted to linearized analysis, which treats the amplitude of 
the unstable long-wavelength magnetic field as perturbatively small. This is adequate for 
calculating the linear growth rate of the instability. In order to understand how instabilities 
affect equilibration, however, we need to know just how big the unstable long- wavelength 
magnetic field grows and what happens next. In particular, does the presence of this non- 
perturbatively large, long- wavelength field cause rapid isotropization and equilibration of the 
system? And what (parametrically) are the time scales involved? These questions cannot 



Whether this will continue to be true in some future bottom-up scenario that incorporates plasma insta- 
bilities from the beginning remains to be seen. 
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be answered by a linear analysis of the instability. In this paper, we will modestly focus on 
the very first question: How big does the instability grow? In answering this question, we 
will also learn important lessons concerning the qualitative nature of what happens to the 
unstable modes of the theory once they have grown large. 

To investigate collective effects, it is both standard and convenient to describe the non- 
equilibrium quark-gluon plasma using kinetic theory in the form of Vlasov equations. We 
split the system into two parts: short wavelength ( "hard" ) momentum excitations which are 
described by a Boltzmann equation, and long- wavelength ( "soft" ) modes which are described 
by classical Maxwell equations. Hard excitations are treated as a collection of particles 
with phase-space density f{p,x,t). For an Abelian theory, the corresponding (coUisionless) 
Boltzmann equation would be 

[dt + v-V^ + g{E + vxB).Vp]f^O, (1.2) 

where v is the velocity associated with p, and E and B represent the fields of the long- 
wavelength modes (those not described by /). The corresponding Maxwell equations for the 
soft modes are 



d.F'^^^f^g / v^f, (1.3) 
Jp 

where there is an implicit sum over particle species on the right-hand side and = {l,v). 
We use the ( — 1-++) metric convention. One can generalize these fully non-linear Vlasov 
equations to non- Abelian plasmas. However, for the discussion in this paper, it is simpler to 
discuss the linearization of the equations in small fiuctuations Sf{p, x.t) of the distribution 
functions about some initial homogeneous distribution fo{p). In this case, 

{dt + v- V,) 5f + g{E + vxB)- Vpfo = 0, (1.4a) 

a.F^"=/ = 5 / v^6f (1.4b) 



p 



for the Abelian theory, assuming that /o carries no net charge or current. The non- Abelian 
generalization which describes the dynamics of long wave-length color electromagnetic fields 
is [20-22] 2 

{Dt + v- D^) 6f + g{E + vx B)- V^/o = 0, (1.5a) 

D^F^^^f^cg I v>^5f, (1.5b) 
Jp 

where /o is colorless and 5f takes values in the adjoint color representation. D is the 
adjoint-representation covariant derivative, and c is a constant that depends on the color 
representation of the hard particles.^ Note that we have linearized in 5f but not in the 



^ Also sec Rcf. [23] for a related formulation. 

^ Specifically, making the implicit sum over hard particle types s explicit, the right-hand side of (1.5b) is 
Y^s^sQ fpV^^fs, where Cg = Vgts = VsCsds/d\. Here, tg is defined by tr(T°Tj') = ts5'^^, where T° are 
the color generators for the hard particle's color representation, Cs is defined by T"T" = Cg, dg is the 

dimension of the color representation, and dp^ is the dimension of the adjoint representation. Vg represents 
the number of non-color degrees of freedom of types s. For instance, for hard gluons in QCD, dg = 8, 
ts = Cs = 3, and Vg — 2 for helicity provided / represents the density of hard gluons per spin state and 
color. 
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strength of the soft gauge field A. This is the usual starting point for kinetic theory discus- 
sions of hard thermal loops (HTLs) [24] in non-Abelian plasmas, whether for isotropic or 
anisotropic /q. 

There are a number of conditions that must be met for coUisionless Vlasov equations to 
give an accurate approximation to the underlying physical situation. The distance scales 
of interest must be (i) large compared to the deBroglie wavelength of the hard particles, 
so that hard particle propagation can be treated classically, and (ii) small compared to the 
mean free path for individual hard particle collisions, so that hard particle collisions may 
be ignored. The soft fields must be large enough that they can be treated as classical fields 
{i.e. the energy in each mode k of interest should be large compared to the frequency ujk of 
that mode). The last assumption will not be a problem for a growing unstable mode, once 
it has grown to significant size. 

In the application to the thermalization of the quark-gluon plasma, the "hard particles" 
are the gluons with momenta of order the saturation scale Qg- The "soft" fields will refer 
to unstable modes of the gauge fields with much smaller momentum. (Note that hard and 
soft are relative terms, and, in much of the literature on saturation, the momentum scale 
Qs would normally be thought of as a "soft" scale.) In order to keep the discussion of this 
paper somewhat general, we will henceforth refer to the hard momentum scale as Phard or 
simply p rather than Qs, and we will refer to the soft momentum scale as fcsoft or simply 
k. We will assume kson <^ Phard- This is the case for the plasma instabilities of the original 
bottom-up scenario at times r 3> l/Qs [19]- 

Though the precise value is not important to our discussion, the constant c in (1.5b) is 
c = 2Ca = 6 for hard gluons if / represents the density of hard gluons per helicity and color 
state. Ca is the quadratic Casimir of the adjoint color representation, with Ca = N for 
SU(7V). 

We can now lay out more precisely the basic question we will discuss in this paper. There 
are two distinct natural scales associated with the question of how big soft fields can grow 
before their effects can no longer be treated perturbatively in Eqs. (1.5). The first has to 
do with the Boltzmann equation (1.5a) and with when soft fields have a large effect on 
the trajectories of the hard particles — that is, with when the assumption 6f -C /o breaks 
down. Consider a soft magnetic field B. The hard particles will follow nearly straight-hne 
{i.e. nearly free) trajectories if the Larmor radius R ~ Pha.id/gB is large compared to the 
wavelength A ~ 1/ fcgoft of the magnetic field. That is, the effects of the soft field on the hard 
particle trajectories are perturbative if 

T) ^soft Phard / „ \ 

-D < . (1.6) 

9 

If one picks a "reasonable" gauge, where the gauge fields A are relatively smooth and do 
not have variation on scales different from /cgoft, then B ~ kA, and this condition can be 
restated as'' 

A<^. (1.7) 
9 



A simple mnemonic for this condition is to think momentarily about the more fundamental gauge theory 
that describes the hard modes, and ask when the soft contribution to the gauge field A in the covariant 
derivative D = d — igA can be treated as a perturbation when applied to the hard modes. When applied 
to hard modes, D phard — gA, and so the condition is gA <C phard- 
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In non-Abelian plasmas, however, we also have to think about the interaction of soft modes 
with each other. A quick way to estimate when these are important is to consider the 
covariant derivative D = d — igA in the non-Abelian Maxwell equation (1.5b). When 
applied to the soft field strength F, we have D ~ ksott — gA. The gA term can only be 
treated as a perturbation when gA <S ksoh, which is 

^«^. (1.8) 
9 

Which of these parametrically different scales determines how large the instabilities grow? 
Is it 

A^^, (1.9) 

giving B ~ kp/g, or 

Ar^^^P^, (1.10) 

9 9 

giving B ~ k'^/g <^ kp/g 7 That is, is the growth stopped by large effects on hard particle 
trajectories, or stopped earlier by non-Abelian interactions between the growing soft modes? 

In the next section, we will briefly review the linearized theory of plasma instabilities. 
In section III, we then evaluate the effective potential for a certain class of unstable con- 
figurations and argue that, when non-linear interactions of the soft modes are considered, 
the instabilities will grow beyond the non-Abelian scale (1.10) by evolving toward Abelian 
configurations. In section IV, we check this assertion in a 1+1 dimensional toy model, 
which we simulate numerically. We find that the growing instabilities indeed evolve into 
gauge configurations living in the maximal Abelian sub-algebra of the gauge theory, which 
is U(l)xU(l) for QCD. In section V, we briefly review what is known about the fate of 
plasma instabilities from studies of traditional plasmas, based on U(l) electromagnetism, 
once they grow large enough to saturate the Abelian bound (1.9). As an aside, in section VI 
we point out that there is a simple lower bound on the exponent n in the parametric relation 
(1.1) for the equilibration time. Finally, we offer our conclusions in section VII. A number 
of topics briefly discussed in the main text are expanded upon in various appendices. 



II. REVIEW OF LINEAR INSTABILITY 

One may formally solve the linearized Boltzmann equation (1.5a) for 6f and then insert 
the solution into Maxwell's equation (1.5b) to obtain an effective equation of motion for the 
soft gauge flelds: 

D.F"" = -c/ [ v^{Dt + v-D^ + e)-\E + v x B) ■ Vp/o, (2.1) 
Jp 

where e is a positive infinitesimal that selects the retarded solution. This equation simplifies 
enormously if one specializes to small gauge fields by linearizing in A. This limit has been 
analyzed by many authors, and the result is similar to that for an Abelian plasma of charged 
particles. Henceforth dropping the subscript on /o, linearization yields 

d^F'"' ~ -cg^ j v^{dt + v-V^ + e)-\E + vxB)- Vp/. (2.2) 
Jp 
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Fourier transforming from {t, x) to K — {cu, k), this can be put in the form 



iK^F^"" ~ -U^^A^ (2.3) 

or equivalently 

Sf'^A^ = 0, (2.4) 

where 

(a;, k) = (-a;' + k^)g^"' - Ri^K" + W{uj, k) (2.5) 

and 

2 /■ df{p) 



U'"'{u;,k)^cg- 



I. 



+ 77 - 

V ■ K — le 



(2.6) 



We will assume that /(p) is parity-invariant in what follows. 

One can now investigate whether there are unstable solutions to (2.4), meaning solutions 

where uo has a positive imaginary part. Ref. [19] uses continuity arguments to show that 
a sufficient condition for the existence of such an instability with a given wavenumber k is 
that the spatial part Sij of S{uj, k) have a negative eigenvalue at zero frequency. That is, 

Si [k'^5ij - kikj + Uij (0, k)] Sj < (2.7) 

for some spatial polarization e (which may be assumed transverse to k without loss of 
generality). This is a generalization of what's known in plasma physics as the Penrose 
criterion. (See also the earlier discussion in Ref. [10].) Roughly speaking, Ref. [19] showed 
that this condition is satisfied for some k when f{p) is anisotropic. Note that a negative 
eigenvalue for Sij{Q, k) requires a negative eigenvalue for the spatial self-energy Ilij{Q, k). 

Ref. [19] also studied in detail a specific case relevant to the early stages of the original 
bottom-up scenario, where f{p) is axially symmetric and extremely anisotropic, with typical 
Pz ^ Px,Py, as depicted schematically in Fig. la. Such a distribution has a flat pancake 
shape in p space. Qualitatively, the resulting set of soft, unstable modes arising from (2.4) 
was found to look like the region of k space depicted in Fig. lb, which has a narrow cigar 
shape aligned along the z axis. In this particular situation, typical unstable modes have fe's 
which are almost parallel to the z axis. Whether such extreme anisotropy will persist in the 
theory of quark-gluon plasma equilibration, once the effect of plasma instabilities are fully 
incorporated into the bottom-up scenario, is not yet known. 

Before moving beyond the linearized effective theory of the soft gauge fields, it will be 
useful to note that the linearized equation of motion (2.4) can be associated with a corre- 
sponding effective action for the linearized theory: 

^(Unearized) ^ f ^_ l _ l ^a^.u^a-^ ^ (2.8) 

J X 

where x stands for (t, x), and J^^u = d^A^ — dj,A^ is the linear piece of the non-Abelian field 
strength F^^. Note that the action is non-local because , given by (2.6), is a non-local 
operator in x. h\ Aq — Q gauge, the action takes the form 

^(Uneaxized) ^ f |^^^|2 _ 1 | ^ x Aj^ - i A.W A,] . (2.9) 
J X 
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(a) (b) 

FIG. 1: (a) Flat pancake shape of the distribution of hard particle momenta p in the early stages 
of the original bottom-up scenario, (b) Associated narrow cigar shape of the set of soft momenta 
k of unstable modes of the soft gauge fields. The scales in the two figures are not the same: 

^soft < Phard- 

Studying the condition (2.7) for instability is equivalent to looking for unstable directions 
of the "potential energy" defined by considering the effective action (2.9) for static (that is, 
a;=0) gauge fields: 

V[A] - I [i|Vx A|2 + iAn,,-(a;=0)^,-] 

= / i Ai{ky - hkj + ni,(0, k)\Aj{k). (2.10) 

Jk 

In (2.10), we have written njj(0, k) instead of njj(0, k). This is because the formula (2.6) 
for n(a;, fe) depends only on the direction k (and not the magnitude) of k when a; = 0. We 
write n(0, k) to emphasize this. In the hmit /c — > 0, the instability condition that S'jj(0, k) 
have a negative eigenvalue is equivalent to 11^^ (0, k) having a negative eigenvalue.^ 

Before continuing, we should mention some limitations of the effective potential. Note 
that, in Aq=Q gauge, the action for a free electromagnetic field J + ||-Bp] looks 
like a kinetic energy term which gives the electric energy, plus a potential energy term 

\\B\'^ which gives the magnetic energy. Generally, the effective potential in = gauge 
represents magnetic potential energy, and electric effects are dynamical effects. As a result, 
one does not see electric effects such as Debye screening looking only at the effective potential 
(2.10). But one can use the effective potential as a way of understanding instabilities that 
manifest as negative eigenvalues of 11^^ (0, k). 



^ We should warn the reader that not all plasma instabilities of the system always correspond to negative 
eigenvalues of n,j(0, fe). See the discussion of "electric" instabilities in Ref. [19]. However, for k pointing 
close enough to the z axis for distributions like Fig. la, the negative eigenvalues of 11^(0, fe) do find all of 
the instabilities [19, 25]. 
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III. THE EFFECTIVE POTENTIAL 



A. Hard particle effects 

In order to investigate whether self- interactions of the soft gauge fields can stop the growth 

of the instabilities, we wish to analyze the effective potential energy without linearizing 
the effective theory in the gauge field A^. Mrowczyhski, Rcbhan, and Strickland [26] have 
derived the effective action which gives the non-linear equation of motion (2.1). Their action 
for anisotropic / is a generalization of the simple form of the HTL effective action derived 
by Braaten and Pisarski [27]^ for isotropic /. It is given by^ 

To get the potential energy, we evaluate this action for static configurations A{x) in Aq = 
gauge, as in the previous section. Even with these restrictions, the formula (3.1) remains a 
rather formal and complicated expression. However, it simplifies tremendously if we restrict 
attention to a certain sub-class of gauge field configurations inspired by Fig. 1. Because 
typical unstable fc's in Fig. lb have A;^; S> generic combinations of these modes will vary 
much more rapidly with z than with x or y. Let us therefore consider the extreme case of 
gauge field configurations that depend only on z: 

A = A{z). (3.2) 

Making use of some observations by Blaizot and lancu [29], we show in Appendix A 1 that, 
for A = A{z), the effective action (3.1) reduces to the simple, local form 

V[A{z)] = / [|i^^i^^ + iA^%(0,e,)A«] 

= [ [lB'^.B^ + lA^U,j{0,e,)A^], (3.3) 

J X 

where 11 is the self-energy (2.6) of the linearized theory, and is the unit vector in the 
z direction. B represents the full non-Abelian magnetic field. This is exactly the same as 
the result (2.10) of the linearized theory applied to configurations A = A{z) except that 
the quadratic | V x term has been replaced by the full, non-Abelian B^, which contains 
cubic and quartic interactions. The term AHA representing the effects of hard particles 
remains quadratic, however, even though we are no longer linearizing the theory in A. 

In field theory, one generally investigates questions of stabihty by studying the effective 
potential for low- momentum modes, taking A; — > 0. If we take this limit in (3.3), we obtain 
the potential energy density 

V = -i/[A,A,-]"[A,A,]« + iA«n,,(0,e,)A^", 

= i g^r'^'^r'^-A^A'^.A^A^. + i A« n,,. (0, e,) A«, (3.4) 



^ The HTL effective action for isotropic systems was originally derived, in a different form, by Taylor and 
Wong [28]. Ref. [26] also discusses the generalization of the Taylor- Wong form to anisotropic systems. 
The difference in overall normalization of (3.1) with Ref. [26], taking c = 2Ca for hard gluons, is due to 
different choices of normalization for /. 
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where / are the usual gauge-group structure constants. We can now investigate the 
topography of this potential in the space of A"'s. We will assume that the hard particle 
distribution functions f{p) are axially symmetric about the z axis and that llij{0,ez) has 
a negative eigenvalue, which is the case for oblate distributions hke Fig. la.^ Using the 
transversality Ki^U.^^'^{K) = of the HTL self-energy (2.6), we can then write 

V - i g'T^J^'^A'^A'^^AfA^^ - i ^^A^A^ + A;A'^), (3.5) 

where 

/i' = -n,.(0, e,) = -UyyiO, e,) > 0. (3.6) 

In the linearized analysis, the unstable modes are those with k < /j,, and the typical unstable 
modes have kgoft ~ (that is, fc < but typicaUy not k <^ /i). 

The potential (3.5) is unbounded below. Specifically, consider any "Abclian" config- 
uration, by which we mean a gauge field where all the components Ai of A commute. 
(For instance, this happens if the gauge field points in a single direction in adjoint color 
space, such as — Ai5'^^.) Then the quartic term in the potential vanishes, leaving 
V = — I ijl'^{A^A'^ + A'^A^), which clearly runs away to — oo as the magnitude of A increases. 
Note that this picture can no longer be trusted when A gets so large that the assumption 
Sf <^ f used to derive the effective action breaks down. That is, the runaway growth of 
Abelian configurations should stop when the size of A is given by the scale Phard/g of (1.9), 
as discussed in the introduction. 

We can get a simple picture of the topography of V for non- Abelian configurations if we 
make some simplifying restrictions on A. We will consider the case where both (i) A lies 
in an SU(2) subgroup of color SU(3) and (ii) A^ = 0. With these assumptions, we show in 
Appendix A 2 that we can always make a combination of spatial and color rotations to put 
A{k — > 0) in the following form: 

Ai = (l)lSiJal + <p2SiySa2- (3.7) 

The potential (3.5) is then 

V(0i, 02) = i g'<l>l4 - \ tJi\<t>l + <t>l). (3.8) 

This potential is depicted in Fig. 2. The Abelian configurations correspond to the 0i axis 
(02 = 0) and the 02 axis (0i = 0). One can see that there exist static non- Abelian solutions, 
indicated by the intersection points of the four straight lines in the figure. At these points, 
the amplitude of the gauge field is A ~ ^/ g. Recalling that unstable modes typically have 
^soft ~ this corresponds to the non- Abelian scale A ~ Koix,/ g discussed in the introduction. 
However, these solutions are unstable to rolling down and subsequently growing in amplitude 
along one of the axis. The picture suggests that, if we start from A near zero, the system 
might possibly at first roll toward one of these configurations with A ~ A^soft/fi', but its 
trajectory would eventually roll away and approximately "Abelianize," growing along either 
the ±01 or ±02 axis until the effective action breaks down at ^4 ~ phard/5'- 

We have been focusing on soft configurations A[z) in the arbitrarily long wavelength limit 
— 0. In Appendix B, we show that the static, unstable, non-Abelian solutions discussed 
above have analogs with finite wavelength as well. These new solutions are also unstable 
and are merely a curiosity: they do not affect the current discussion. 



For a qualitative discussion, see Ref. [19]. For calculations in various cases, see Refs. [13, 25]. 
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FIG. 2: Two equivalent depictions of the potential V{(f)i,(j)2) of (3.8). The (pi and 02 axis are in 
units of fi/g, and the values of V are in units of l^^/g^- There is a local maximum with V = 
at the origin. The four straight lines in the plots correspond to the equipotential V = —^'^/2g^, 
and the four intersection points of those lines are saddle points corresponding to static, unstable, 
non-Abelian configurations. 

B. Soft mode effects 

So far, we have constructed the effective action and the effective potential by integrating 
out the effects of the hard particles. We have not yet, however, completely considered the 
effect that various soft field excitations with different fc's can have on the instability. In 
order to illuminate the remaining issues, let us consider a very simple, warm-up toy model 
for the effective Lagrangian for the soft fields, which is a scalar field theory with the potential 
V(0i,02) of (3.8): 

£ = -\{dct>if - \{dct>2f - V(0i, 02). (3.9) 

As an example of a concern one might have about our previous discussion of instability, 
imagine that the soft sector were at some finite temperature T ^ /z. What happens to our 
discussion when we account for the effects of interactions with modes with momenta of order 
T? It is well known that such interactions give a contribution of order g^T"^ to the effective 
scalar mass-squared m? . This effective Am? is generated though diagrams such as Fig. 3a, 
which physically represent the forward scattering of particles off of the thermal bath, as in 
Fig. 3b. As a result, 

-fi^ ^ -fi^ + 0{g^T^) (3.10) 

in the effective potential (3.8). If T were large enough, this would stabihze the potential 
near the origin and prevent runaway growth of the fields.^ This is the standard picture of 
high-temperature symmetry restoration in scalar field theories [30-32]. 

Can anything similar happen in the problem of interest, where we consider an initial 
state that has very little energy in soft modes, but these modes subsequently grow because 
of the instability? In our toy scalar model, the answer is no. For the scalar theory, we give 
simple qualitative arguments (in any number of dimensions) in Appendix C. But a more 



^ Such a thermal mass effect would in fact only render the origin me<a-stable in this scenario. The effective 
mass of small fluctuations behaves as gcf> for large (j> and exceeds T for ^ T/g. Thermal effects are then 
suppressed by exp(— g^/T), and so the stabilizing thermal contribution to the mass would disappear at 
large enough (p. 
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(a) (b) 



FIG. 3: (a) 1-loop mass correction in finite-temperature diagrammatic perturbation theory, (b) 
Forward scattering off of a particle in the plasma. 

straightforward way to verify our assertion is by explicit numerical simulation, which we 
discuss in the next section for 1+1 dimensions with a more sophisticated toy model than 
the scalar theory (3.9). 

IV. NUMERICAL ANALYSIS OF ABELIANIZATION IN A TOY MODEL 

There are several uses for performing a numerical simulation of a representative toy model 
of the physics we have been discussing. One is simply to verify our claim in Sec. Ill B that 
soft excitations produced during instability growth cannot stabilize the system from growth 
beyond the non-Abehan scale (1.10). Another is to better understand whether Abehaniza- 
tion is a global phenomena. We have suggested that the system, in searching to minimize its 
potential energy as the fields grow large, seeks configurations where commutators [Ai, Aj] 
are small, since these would otherwise contribute substantially to the non-Abelian magnetic 
field energy B'^. But is this a local statement or a global statement? Can one necessarily 
find a gauge where all Ai{x) commute with all Aj{y) for well-separated x and yl That 
is, over what distance scale can one say that the non-Abelian gauge configurations are well 
approximated by purely Abelian gauge configurations, and does this distance scale grow 
with time? In this section, we will address these questions by numerical simulations of the 
real-time evolution of a 1-1-1 dimensional gauge theory. 

A. dimensional gauge theory toy model 

Consider the full effective theory (3.1) of soft modes but restrict attention to gauge fields 
of the form 

A^ = A,{z,t). (4.1) 

This is similar to the restriction A = A{z) of the last section, except that we now allow 
for time dependence. In addition, we will ignore time dependence in the second term of 
(3.1) so that it becomes the simple local ^A^llij(0, ez)Aj of (3.5). This last approximation 
is not in general justified, which is why this will only be a toy model calculation, and we 
will discuss some of its deficiencies in Sec. IV E. We hope it will qualitatively capture the 
physics of interest. In 1-1-1 dimensional language, Aq and A^ are the gauge fields, and A^ 
and Ay behave as adjoint-representation scalars. To emphasize this distinction, we will refer 
to Arc and Ay as (px and (py in what follows. The toy model Lagrangian is then 

a=x,y 
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where V{(l)x,4>y) is the potential of (3.5). By writing = T"0", where are fundamental 
representation color generators with the usual normalization tr(T"r'') = we can also 
write this in the form 

jC = tT{El) + tr[(L>o0,)2 - {DJ,f] + tr[(L>o0,)^ - {D,<Pyf] - V(0,, 0,), (4.3) 

V(0., <Py) = V tr(0^ + 05) + g' tr ((^[0,, 0,])') . (4.4) 
In = gauge, the equations of motion are 

^1 = {Dl + i^^)cl), - /[02, [02, 0i]], (4.5a) 

02 = (^,' + /i')02 - [01, 02]], (4.5b) 

A = i5r[0^,i:)20^] +ig'[0y,D^0j,], (4.5c) 

with 1^20 = (9^0 — ig[Az, 0]. Gauss's Law, which is a constraint equation preserved by the 
equations of motion, is 

D^Az = ig[(l)^, 0a;] + ig[(l)y, (j)y]. (4.6) 

Classical evolution does not depend in any essential way on the value of g. One can 
remove g hj a simple rescaling of fields, A^/g (including A^ = (px and Ay = (j)y). This 

is equivalent to setting g = 1 in the equations (4.5). Recall that the typical unstable modes 
have momenta of order fx. It is therefore also natural to work in units where = 1, so that 
the condition (1.10) for when non-Abelian interactions between growing unstable modes first 
becomes important is simply ~ 1. However, we will retain factors of g and n in what 
follows simply so that our notation and discussion is consistent with the conventions used 
earlier in this paper. Readers are encouraged to ignore these factors if so inclined, setting 
g — 1 and /j, ~ 1 in what follows. 

Our goal will be to evolve the system (4.5) classically, starting from small, random initial 
conditions for the fields. Note that Aq=0 gauge is preserved by time-independent gauge 
transformations. In particular, in infinite volume, ^'^ one can always find a time-independent 
gauge transformation to put the fields in axial gauge = at a particular time ti. We 
will use this freedom to choose our initial condition at t = to have — 0. We want 
tiny initial fluctuations in our soft flclds to provide seeds for the soft-mode instabilities we 
wish to simulate. For A^ = (px and Ay = (t)y, we choose their initial values to be zero and 
their initial time derivatives to be Gaussian random white noise with a very small amplitude 
A.^^ Gauss's Law (4.6) then determines that — —A^ must be z independent at t = 0. 
Since we are not interested in background electric fields, we take to be initially zero. In 
summary, our initial conditions (written in continuum notation) are 

A,{z, 0) = 0,(z, 0) = (t>y{z, 0) = 0, E,{z, 0) = 0, (4.7a) 
(0«(z,O)05(/,O)> ^^5'^'5^p5{z-z'). (4.7b) 



We will be simulating finite volumes, where there can be non-trivial spatial Polyakov loops, but we use 
the infinite- volume case to inspire our choice of initial conditions. 

Some readers may worry that using a white noise distribution may populate UV lattice modes with enough 
energy to seriously distort, through interactions, the physics of the soft sector in the continuum limit. 
This is not a problem and is discussed in Appendix D. 
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In our simulations, wc work on a 1 -dimensional periodic spatial lattice with lattice 
spacing a = 0.05 /i^^ and length L = 250 fi^^ = 5000 a. To evolve the system in time, 
we use a leapfrog time algorithm with time step e = 0.001 /x^^. Our algorithm very 
closely follows that of Refs. [33, 34].^^ Our discrete-time lattice evolution equations, 
along with the discretized version of our initial conditions, are given explicitly in Ap- 
pendix D. We choose A, which parametrizes the size of the initial fluctuations of 0, to 
be A = 0.01 {anY^^fx^/g ~ 0.002236 /x^/^f. We evolve a single, representative, random initial 
field from the ensemble defined by (4.7b). 

In addition to simulating the model for SU(3) gauge theory, we have also simulated it for 
SU(2) gauge theory. The color structure of SU(2) is slightly easier to discuss pedagogically 
because its maximal Abelian subgroup is U(l). If two adjoint fields commute in SU(2), they 
must be pointing in the same color direction. SU(3) is more complicated. For this reason, 
we start by discussing SU(2). 



B. Simulation Results: SU(2) gauge theory 

Fig. 4 shows how the amplitude of the fields grows with time. Specifically, it shows the 
growth of the rms average 



101: 



/ f + - [/ ^2tr(0^ + 0^) 



1/2 



(4.8) 



The figure also shows the volume average of the relative size of the local commutator 
of the fields = (j)^ and Ay = (py, as measured by 



C 



dz 



C{z), 



C = 



tr 



(^{i[(j)x,(f)y]y 



1/2 



tr(02 + 



(4.9) 



As one can see, there is a stage of the instability growth, around ~ where the 
original, random, non- Abelian fluctuations suddenly change character and (px and (j)y become, 
at least locally, aligned in the same color direction. We will analyze soon how global is this 
alignment. 

Before continuing, we should note that quantities such as tr(A^ + Ay) and tr([A2., A^,]^) 
are gauge- invariant under 1-1-1 dimensional gauge transformations, but they would not be 
be gauge-invariant in the original 3-1-1 dimensional theory under general 3-1-1 dimensional 
gauge transformations. This makes the construction of gauge-invariant observables easier in 
the 1-1-1 dimensional theory than in 3-1-1 dimensions. In particular, we can directly probe 
statements about the relative size of the commutator [A^, Ay] without having to either (i) 
rely on the vague, qualitative statements like "in a reasonable gauge" which we used earlier 
in this paper, or (ii) construct more indirect observables, like commutators of magnetic 
fields. We shall take advantage of this feature of the 1+1 dimensional theory. 

We have shown what happens to the relative size of [0^, 0^]. Readers may wonder about 
the other possible commutators of the basic fields of our model, [^42,0^.] and [^4^,02^]. The 
value of [Az, 0q] at a particular moment is not physical because we can make it vanish 
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See also the overview in Ref. [35] and the closely related algorithms discussed in Refs. [36, 37]. 
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FIG. 4: The average amplitude |0|rms (solid) and the relative size C of commutators (dashed) as 
a function of time for gauge group SU(2). 

everywhere by gauge transforming so that /l^ = at that moment. We can think of fixing 
axial gauge at a particular time t as the 1+1 dimensional analogue of a "smooth" gauge 
choice for that time slice. 

Snapshots of the field configurations are shown at five different times in Fig. 5. The first 
row shows the field strength 

$ = {rAl + Wyf'^ = [2 tr(0^ + ct>l)\ (4.10) 

as a function of z. The second row shows the relative size C of commutators, as defined in 
(4.9), again as a function of z. Note that the scales used in the graphs change with time. 

We have learned that (px and (j)y point in the same color direction at late times. We would 
now like to address whether this direction "changes" with z. That is, is it possible to find 
a gauge where the fields all live in an Abelian sub-algebra of SU(2), e.g. all pointing in the 

direction: 

Al{z, t) = A,{z, t) C(z, t) = <Po.{z, t) (4.11) 

A gauge-invariant way to compare the color directions of (j)a{z) for different z is to parallel 
transport all of the (f)a{z) to some reference point Zq. We will pick Zq = 0. Equivalcntly, one 
can simply gauge transform to A^=0 gauge on the particular time slice of interest, in which 
case parallel transport is trivial and the color directions at different z can be compared 
directly. 

In SU(2), the color direction of an adjoint field can be represented as a unit 3- vector 
(0^0^03)/(0a^a^l/2_ therefore lives on a 2-sphere S"^. However, the fields of an Abelian 
configuration (4.11) can take both positive and negative values, corresponding, for example, 
to color directions (0, 0, ±1). So, to test how much a given configuration deviates from being 
Abelian, we should not differentiate between a given color direction and its negative. The 
relevant color space is therefore with antipodal points identified (6*^/^2). The third row 
of Fig. 5 shows the (parallel transported^^) color directions swept out as z varies from 



In a periodic space, there are two different ways to parallel transport to ^; = 0: to the left or to the right. 



14 



nt=l 



iit=5 



Iit=10 




0.002 




50 100 150 200 250 
\i.Z 



0.4 
0.2 



50 xOO 150 200 250 






50 100 150 200 250 50 100 150 200 250 




I * 0.2 



50 100 150 200 250 50 100 150 200 250 





Ht=12 




Ht=14 



1200 
1000 
800 
600 
400 
200 




50 100 150 200 250 



50 100 150 200 250 



0.175 
0.15 
0.125 

C 0-1 

0.075 
0.05 
0.025 










007 













006 













005 













004 













003 













002 






.. J .1 1 L ll 







001 


. 1 i .1 1 1 J 












,1 1 



50 100 150 200 250 



50 100 150 200 250 



FIG. 5: Snapshots in time of the evolution of SU(2) field configurations. Each column corresponds 
to a different time, labeled at the top: fit = 1, 5, 10, 12, 14. For each time, the first row shows the 
total field strength (4.10), in units of /i/g, as a function of fj,z. The second row shows the relative 
size C of commutators as a function of ^z. The third row shows the the color directions swept out 
in axial gauge as one varies z from to L, as described in the text. The circles mark the points on 
these curves where /iz is a multiple of 50, and squares indicate the endpoints z = and z = L. 
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FIG. 6: Projection of the space S'^ of SU(2) adjoint color directions into the plane, with the 
identification of antipodal points. 



FIG. 7: The angle in the xy plane of the vector (f) = {(f)x,(t)y), as a function of fiz at /xt = 14, for 
SU(2) simulations (taking (p projected onto the color direction indicated in the last plot of Fig. 5). 

to L. The color direction^^ is represented in the following way. Consider a color direction 
represented as a point on the unit sphere, as shown in Fig. 6. If it is in the lower hemisphere. 



Our convention in these plots is to always parallel transport to the left. Equivalently, there may be an 
obstruction to transforming to = gauge everywhere in a periodic space. Our convention here is 
to transform to = everywhere except on the link that connects across the periodic boundary that 

identifies z = and z = L. 

Readers may wonder which color directions we have plotted, since there are two different fields, 4>x and 
For an Abelian field, = {(px,4'y) points in a single, well-defined direction £{z,t) in the xy plane 
at each point. For general non- Abelian fields, we simply define e as the unit vector which maximizes 
(e • (f>)"'{e ■ 0)". The plots in the third row of Fig. 5 are of the color direction of e • 0. At late times, this 
aligns with the color directions of cp^ and (j)y except at those isolated points where e is perpendicular to x 
or y (where the very tiny non- Abelian components of (jj^ or (py will dominate over the Abelian one, which 
is simply an accident of how one chooses the x and y axis). 








50 100 150 200 250 
|i z 
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use the identification of antipodal points to replace it by a point in the upper hemisphere. 
Then project the points in the upper hemisphere down to the plane passing through the 
equator. The plots in Fig. 5 show that plane. The dashed circles represent the equator, for 
which antipodal points are identified. 

One can see from the color direction plots in Fig. 5 that the color directions globally align 
as t gets large. At fit = 14, the configuration is to very good approximation Abelian over 
the entire simulation volume. However, other aspects of the field remain uncorrelated. Fig. 
7 shows the angle corresponding to the direction in the xy plane of the essentially Abelian 
</> = (0a;, (f)y) at lit — 14, as a function of z}^ So, the SU(2) configurations at late times are 
homogeneous in color and look like a random superposition of unstable Abelian modes. 

C. Abelianization Length: results for SU(2) 

The adjoint color space of SU(3) is more complicated than that of SU(2). When we study 
SU(3) in the next section, we will not be able to make simple, visual plots of color direction 
as in Fig. 5. Also, the maximal Abelian subgroup of SU(2) is U(l), but the maximal Abelian 
subgroup of SU(3) is U(l)xU(l). In SU(3), colors need not point in a single color direction 
in order to be Abelian. It will therefore be useful to replace our pictorial analysis by an 
appropriate correlation length that measures over what distance scales fields commute with 
each other. Consider the following correlation, which measures whether fields commute over 
a given distance ^: 

dA dz ['^"(^ + 0' ^(^ + ^' ^) ' ['^«(^ + 0, U{z + e, z) (t)p{z)\ 



I 

Jo 



(4.12) 



'2CaJo L tr{(l)^(z + ^)(l)a(z + ^)}tr{Mz)Mz)} 

where U{z', z) represents adjoint-representation parallel transport from z to z' . Specifically, 

U{z\ z) M^) = U{z', z) M^) [U{z', z)] \ (4.13) 
where U{z',z) is the fundamental-representation transporter 



U{z',z) =Pexp 



tg f dz"A,{z") 

J z 



(4.14) 



and V indicates path ordering (with z' on the left and z on the right). Alternatively, one 
may transform to A^—^ gauge and dispense with the parallel transporters.^^ In (4.12), d^ 
and Ca are the dimension and quadratic Casimir of the adjoint color representation. 

If commutators vanish locally, as we have seen they do for late times, then Xa(0) = 0. 
The correlation xa has been normalized so that if the colors are completely uncorrelated 



More precisely, we plot the direction e defined in footnote 14. 

There is a subtlety to this on the lattice when z < L < z + ^ because one cannot fix Az=0 gauge on one 
link, which we choose as the link across the periodic boundary. One can either incorporate this single link 
in the parallel transport, or else extend the axial gauge transformation to the periodic copies of the fields 
beyond z = L. 
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FIG. 8: The solid lines and circles show the Abelianization correlation length as a function of 
time, measured every 5t = 0.1 Individual points, rather than a continuous line, are shown 
simply because jumps around quickly once it becomes a sizable fraction of the lattice. The 
dashed line shows the ordinary correlation length (^(^. The topmost axis of the outer graph is 
L/2 and represents correlation lengths which exceed the size of our lattice. The inset provides an 
expanded view of the behavior for /it < 10. 



over a distance ^, then xiO — 1- We define the "Abelianization correlation length" ^a as 
the smallest distance for which Xa(0 ^ A plot of .^a vs. time for SU(2) is shown in 
Fig. 8. When ,^a exceeds L/2 (indicated by the top of the plot), then the correlation does 
not drop below 0.5 anywhere on the lattice. One can see that this correlation length begins 
to grow rapidly at roughly the same time the relative size C of commutators begins to drop 
in Fig. 4. 

For comparison, we also plot in Fig. 8 a correlation length defined in terms of the full 
correlation (not just color) between the parallel-transported fields, defined by 



f t r{Mz + OUiz + ^,z) 0„ (z) } 



This correlation is normalized so that x(0) = 1, and xiO vanishes if the fields are uncorre- 
lated over distance ^. We define a correlation length by when this correlation first drops 
below 0.5. Note that this correlation length does not grow enormously like the Abelianization 
length ,^A does. 



D. Simulation Results: SU(3) gauge theory 

We are now ready to discuss simulation results for the QCD version of our 1+1 dimen- 
sional model. Simulation results for |0|rms, C, and for SU(3) are shown as a function 
of time in Figs. 9 and 10. The results are qualitatively similar to the SU(2) case, and we 
can draw the same conclusion from the rapid growth of ^a: the configurations Abelianize as 
the field strength grows large compared to the natural non-Abelian scale A fi/g. 
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FIG. 9: As Fig. 4 but for SU(3) gauge theory. 




FIG. 10: As Fig. 8 but for SU(3) gauge theory. 



Plots of $ and C vs. z are qualitatively similar to those of the SU(2) simulations shown 
in the first two rows of Fig. 5, and we will refrain from displaying them. 

The Abelianization in SU(3) is slightly different than in SU(2) because the fields do 
not have to point in a single color direction. The maximal Abelian subgroup of SU(2) is 
U(l) xU(l) and, after an appropriate gauge transformation, can be considered to be spanned 
by the diagonal generators 




(4.16) 



in the standard Gell-Mann representation. The analog of (4.11) for Abelian configurations 
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FIG. 11: The parameters rjx (solid) and rjy (dashed) of (4.18) vs. z for fit = 13, showing the spatial 
variation of color within the maximal Abelian subgroup of SU(3). 

is that they can be put in the form 

A:{z,t) = Al{z,t)6-' + Al{z,t)6^', C{z,t) = <Pl{z,t)6'^-' + <Pl{z,t)6^'. (4.17) 

To investigate whether the late-time SU(3) configurations vary in color with z within the 
maximal Abelian subgroup, we would like to plot the relative and components as 
a function of z. and are somewhat artificial choices, however, because they break 
the permutation symmetry of the U(l)xU(l) subgroup, treating the first two colors of the 
fundamental representation differently from the third. In Fig. 11, we plot instead the gauge- 
invariant measures 

, , 3-\/6 det (j)aiz) , . , . 

ria(z) = -TTT no sum on a) 4.18 

vs. z. These parameters vary in magnitude |?7| between (i) zero, for color matrices with 
eigenvalues proportional to (1,-1,0), like and its permutations, and (ii) 1, for color 
matrices with eigenvalues proportional to (1, 1, —2), like and its permutations. A plot of 
rja{z) vs. z is given for time nt = 13. As one can see, the color direction within the maximal 
Abelian subgroup remains uncorrelated at this late time, even though the Abelianization 
correlation length is larger than the size of the lattice. 



E. Differences between the toy model and the full effective theory 

Before leaving our discussion of simulation results, we should mention one of the qualita- 
tive differences between our toy model (4.5) and the real theory of plasma instabilities. In 
constructing the toy model, we ignored time dependence in the second term of the effective 
action (3.1). Consider a linear analysis of instabilities around zero field. In our toy model, 
the growth rate of each mode with k < fi is given by 

7toy(A;) = (4.19) 
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FIG. 12: Qualitative behavior of linear growth rates of instabilities vs. momentum k for (long- 
dashed) the toy model, (solid) actual QCD with a generic axi-symmetric hard particle distribution 
/(p), and (short-dashed) actual QCD with a pancake-shaped hard particle distribution. 

which is easily derived by linearizing the equations of motion (4.5). Qualitatively, this 
growth rate has the form of the long-dashed line in Fig. 12 and is 7toy = /x at A; = 0. In the 
actual theory, where time dependence is not ignored, 7(A;) vanishes at A; = [10, 13, 19, 25]. 
For generic hard particle distributions, 7(A;) then has the qualitative form of the solid line 
in Fig. 12. The slow response of the system for small k is related to Lenz's Law: changing 
magnetic fields create electric fields, which induce currents to oppose the change. For the 
extremely anisotropic case of pancake distributions, such as shown in Fig. 1, instability 
growth is relatively slow compared to ii for all unstable k [19].^'' This is depicted by the 
short-dashed line in Fig. 12. We find it plausible that slow down of the growth rate will not 
affect the qualitative conclusions we have drawn: that energetics favors the development of 
Abelian configurations. However, this is something that should be checked in the future by 
simulations of the full theory. In particular, one might wonder whether the suppression of 
the growth of small k modes in the full theory might affect the growth of the color correlation 
length to scales large compared to ji"^. 



V. FATE OF ABELIAN INSTABILITIES 

If QCD plasma instabilities indeed Abelianize, then we may be able to learn qualitative 
lessons about their behavior from studies of traditional electromagnetic plasmas. Indeed, 
for SU(2) gauge theory, we have seen in our toy model that plasma instabilities grow into 
configurations of U(l) gauge theory. For QCD, the situation is a little more complicated: 
we get two copies of "electromagnetism" as U(l)xU(l). Still, we may hope to gain some 
qualitative insight from what's known about U(l) gauge theories. In this section, we will 
review a few relevant results from the traditional plasma literature. 

Throughout this paper, we have focused on configurations which are approximately in- 



This follows from the results of Ref. [19] with the realization that <C gmax in that reference when the 
hard particle distribution is extremely pancake shaped. The ^max of Ref. [19] is the ^, of this paper. 
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FIG. 13: (a) Qualitative depiction of a non-linear magnetic BGK wave. The lines represent mag- 
netic field lines. The points and crosses represent the net charged particle current, pointing cither 
out of or into the page, (b) A magnetic tear or reconnection instability for the central region of 
figure (a). 



dependent of x and y, which arise naturally, for example, when the distribution of hard 
particles is pancakc-shapcd. One might wonder whether the dynamics continues to main- 
tain the approximate independence of the configurations on x and y once the instabilities 
grow large enough that their effects on the hard particles become non-perturbative. Could 
the system settle down into oscillations about some non-linear z-dependent magnetic config- 
urations, such as shown in Fig. 13a? Researchers have found examples of time- independent, 
non-perturbative wave solutions of the full (Abelian) non-linear collisionless Vlasov equa- 
tions, such as the magnetic Bernstein-Greene-Kruskal (BGK) wave solutions discovered by 
Davidson et al. [38, 39]. However, one might guess that such solutions are inflicted with yet 
another type of plasma instability, known in different contexts as magnetic tear or recon- 
nection instabilities."*^^ 

Consider the magnetic field configuration depicted in Fig. 13a. By Ampere's Law, a 
stable, time-independent solution of this form requires currents in the x direction, also 
shown in the figure. Visualize the currents as being carried by wires. Two parallel wires 
with current in the same direction attract each other through magnetic interaction. There 
can then be an instability for the wires to dimerize — that is, clump together, as shown in 



Different nomenclature is used by different people. Usually, these instabilities are studied in the magneto- 
hydrodynaniic limit (MHD), which is the opposite limit of the collisionless plasmas studied in this paper, 
as MHD applies to physics on scales large compared to the (transport) mean free path. However, there 
are analogous processes in the collisionless limit. See, for instance, the analysis in Sec. 6.2.2 of Ref. [40]. 
When reading the plasma literature, however, it is important to keep in mind the difference between the 
MHD and collisionless limits since the physics can be different. 

The regions of current in such situations are referred to in the plasma literature as "current sheets." 
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Fig. 13b. The magnetic fields then change correspondingly, as indicated in the figure. In 
the context of collisionless plasma theory, an analytic analysis of this instability in a similar 
situation may be found in Ref. [40].^° It has been observed in numerical simulations of 
non-relativistic plasmas by Calif ano et al. [42], which we will review in a moment. The 
interesting qualitative feature of the tear/reconnection instability is that it breaks the xy 
translation invariance! The non-perturbative dynamics of fully grown plasma instabilities 
therefore has the potential to isotropize the system on relatively short time scales. 

Califano et al. study the case of a hard particle distribution corresponding to two similar, 
uniform, counter-streaming beams of electrons which have a tiny thermal spread of velocities 
around the beam velocities, as depicted in Fig. 14a. Additionally, rather than taking 
random initial conditions for the seed magnetic fields, they take initial conditions of the 
form 

B — (small) X sm.{kQz) By + (much smaller random noise). (5.1) 

The first term seeds a single unstable mode that is xy translation invariant and grows into 
something like the non-hnear magnetic BGK mode discussed earlier. The second term is even 
smaller, and seeds the subsequent tear/reconnection instability of the magnetic BGK mode. 
They made the two terms different sizes because they wanted to clearly see a magnetic BGK 
mode and its subsequent demise, rather than having it all jumbled together. A qualitative 
depiction of the subsequent evolution, pieced together from Refs. [38, 42], is shown in Fig. 
14b. Tx, Ty, and Tz are measures of the total kinetic energies along the x, y, and z axis 
respectively, defined in terms of mv'^, mVy, and mv'^. Initially, almost all of the electron 
kinetic energy is in the x direction. Due to the Weibel instability, the first term of (5.1) then 
seeds the growth of a single, dominant unstable mode. The B field of that mode points in 
the y direction and so can defiect the particles in the xz plane but does not affect Vy. The 
growth of this instability eventually saturates in a nonlinear magnetic BGK-like state. Then, 
at late times, the second term of (5.1) grows due to the tear/reconnection instability. This 
instability drives the system toward isotropization of T^,, Ty, and T^. Unfortunately, Califano 
et a/.'s simulations end before one can see for sure whether there is complete isotropization. 

One may wonder whether plasma instabilities always isotropize a collisionless plasma. 
Refs. [43-46] numerically study parity-asymmetric initial conditions which do not isotropize. 
In particular, Honda et al. [43] study the initial condition of a homogeneous beam of fast 
electrons in a plasma of slow electrons (and slower ions). As time evolves, they find that the 
Weibel instability causes clumping of the beam in the transverse direction. The resulting 
current filaments then attract each other magnetically (again just as two parallel wires 
carrying current in the same direction), and the filaments begin to merge, creating larger 



Bergcr and Davidson [39] perform a stability analysis for a class of magnetic BGK solutions. However, 
they only study stability within the subspace of xy translation invariant configurations. As a result, their 
analysis would not find a tear/reconnection instability. For the same reason, the tear/reconnection insta- 
bility cannot appear in the original simulations of Davidson et al. [38] , nor in more recent one-dimensional 
simulations by Yang, Arons, and Langdon [41] of the Weibel instability in relativistic electromagnetic 
plasmas. 

Here and throughout, we have translated the coordinate labels used by Califano et al. [42] from (a;, y, z) to 
(a;, z, y), to make the dominant term in (5.1) depend only on z, as has been our convention in this paper. 

Their simulations arc 2D+3V, meaning they treat 5 dimensions (2 space and 3 velocity) of 6-dimensional 
phase space. That is, they take configurations to be homogeneous in one spatial direction. 
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(a) (b) 

FIG. 14: (a) A two-stream initial condition for /(p). (b) A qualitative picture of the resulting 
evolution of the initial magnetic seed field (5.1). The right half of figure (b) is a smoothed, 
qualitative sketch of the numerical results of Ref . [42] . The left-hand side was not shown specifically 
in Ref. [42] but is pieced together from statements in the text of that reference and from earlier 
simulations of the Weibel instability in Ref. [38]. 

and larger filaments. Finally, there is only one filament left in their simulation volume, and 
the plasma stabilizes in this configuration. This final configuration is clearly anisotropic and 
corresponds to an equilibrium solution to the Vlasov equations originally found by Bennett 
in 1934 [47], which we review for the ultra-relativistic case (where it is simpler) in Appendix 
E. 

Because of the parity-noninvariance of the initial state, the simulations of Honda et al. 
could end up with a parity-noninvariant final state, described by the Bennett self-pinching 
filament. It is possible, in contrast, that generic parity- muarian^ initial conditions may lead 
to isotropization through coUisionless processes, as suggested by the results of Cahfano et 
al. 

VI. LOWER BOUND FOR THE COMPLETE THERMALIZATION TIME 

In the introduction, we pointed out the embarrassing fact that theory has so far been 
unable to determine even the power n in the parametric relation Teq ~ a~'^Q~^ for weak 
coupling, where Tcq is the local thermalization time. In this section, we point out that there 
is a simple lower bound on n. Consider the last part of thermalization, where the system is 
approximately but not fully equilibrated. At this point, we can use near-equilibrium results 
for the time scales associated with approach to equilibrium. Roughly, the characteristic 
time scale is then simply the equilibrium time scale for (i) large-angle deflections of a given 
particle's direction, and (ii) number- changing processes such as hard gluon bremsstrahlung 
and qq creation or annihilation. Both of these processes have been considered in detail in 
the literature analyzing near-equilibrium transport processes, such as the calculation of the 
quark-gluon plasma shear viscosity [48-50]. Parametrically, the time scale for both these 
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types of processes is 



(6.1) 

ail 



up to a factor of ln(l/Q;s), which we will not bother to keep track of. Full local equilibration 
cannot happen any faster than this, and so 

r^^^- (6.2) 

If the system has fully equilibrated at time Teq during the 1-dimensional expansion, then the 
energy density at that time. 



must be T^, so that 



T - 7-TS%TT74 • (6-4) 



Combining (6.2) and (6.4), we obtain the lower bound 

QsTeq > af'\ (6.5) 

Note that the result Q^Teq ~ Os"^^^^ of the original bottom-up scenario [7] satisfies this 
bound. In that case, thermalization happened slower than the bound because the energy of 
the system got hung up for a time in hard modes, before it could equilibrate by cascading 
down into softer modes. The effects of plasma instabilities might speed up equilibration, 
and it is conceivable that the real answer could saturate the lower bound. Even if it does, 
one would still like to know all the various time scales associated with the different stages 
of equilibration, and what happens at each stage. The time scale for complete local equi- 
libration might not be the most relevant time scale for understanding the physics of heavy 
ion collisions. 



VII. CONCLUSIONS 

In this paper, we have examined the effect of non-Abelian self-interactions on growing 
plasma instabilities in an anisotropic non-Abelian plasma. By looking at the effective po- 
tential for z-dependent magnetic fields in the presence of an anisotropic distribution f{p) of 
hard particles, we have given suggestive arguments that non-Abelian interactions drive the 
growing instabilities to become Abelian once they grow large enough. Numerical simulations 
of a toy model 1+1 dimensional gauge theory also show this behavior. We conjecture that 
(i) the plasma instabilities of SU(2) gauge theory grow into those of traditional (relativistic) 
U(l) gauge theory, and (ii) the plasma instabilities of QCD grow into those of U(l)xU(l) 
gauge theory, which is like traditional plasma theory but with two copies of electromag- 
netism. In the introduction, we asked what sets the scale for how large plasma instabilities 
grow. Our proposed answer is the Abelian one: the scale A ~ Psoft/g for hard particles 
to be affected non-perturbatively, rather than the scale A ~ kgott/g for soft modes to have 
significant non-Abelian interactions. 
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Clearly, these conclusions need to be verified in 3+1 dimensional models, and with the 
full HTL effective action (3.1). In order to better understand the development of plasma 
instabilities once they are grown, full simulations of the full non-linear Vlasov equations are 
needed. Further simulations of traditional U(l) theories would provide useful information. 
Simulations of U(l)xU(l) theories would be interesting as well. Of course, full simulations 
of non-linear SU(3) Vlasov equations would be even better. It is also possible that non- 
equilibrium QCD plasma physics could be studied in simulations of pure classical gauge 
theories on the lattice, in situations where there is a momentum scale separation so that 
excitation of some modes can be thought of as hard "particles" and others as soft fields. 
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APPENDIX A: THE EFFECTIVE POTENTIAL FOR A ^ A{z) 
1. For A = A{z) 

In this appendix, we will derive the result (3.3) for the effective potential when A = A{z). 
We will begin, however, with a more general discussion. We find it advantageous to use a 
slightly different form of the effective action (3.1): 

Ses^- [ \F%,F->"' -eg'' I [ ^ W\W"^, (Ala) 

J X J X J p P 

where 

Wa - Wa{x, v) = ^ F^^{x). (Alb) 

This form can be obtained from (3.1) by using the anti-symmetry of the operator v ■ D in 
x/color space to move one factor of {v-D)~^ from one F to the other. Here and throughout 
this appendix, we will not bother keeping track of the e prescriptions for retarded behavior. 
We note in passing that, in the isotropic case, lancu [51] has discussed a Hamiltonian 
formulation which is similar to (Al). 

Evaluating this action for static configurations A{x) m. Aq — {) gauge, we obtain 

Kff = / \F^F^ + ^V, (A2a) 

J X 



In more detail, the adjoint-rcprcscntation operator D is a real anti-symmetric operator in x/color space. 
That means that the inverse {v ■ D)^^ of f • _D is anti-symmetric as well. Now think of operating on some 
state \s) in configuration/color space, which represents a real function s°'{x) with a single adjoint color 
index. Then the anti-symmetry can be used to rewrite {s'\{v D)-'^\s) = -[{v D)-^\s')] {vD)-^\s). 
Taking s' and s of the form and v^'F^J^ = -v"F^"^, one can then obtain (Al) from (3.1). 
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where 

AV = cg' f W^W^, (A2b) 

and now 

Wk = Wk{x, v) = F,k{x). (A2c) 

Now we specialize to A = A{z). One can then make use of an identity noted by Blaizot 
and lancu [29], which is that 

v.d(a''- ^ (A3) 

\ n ■ V J 

HA — A{n • x) is only a function ofn- x, for some constant 4-vector n. This identity can be 
checked by explicitly expanding all the terms on both sides. Now apply {v • D)~^ to both 
sides of (A3) to get 

n°'n ■ A 

W'^^A"- . (A4) 

n • V 

In our case of A = A{z) with = 0, the identity is 

Wk^Ak-Sk,"^. (A5) 

Vz 

Substituting into (A2b), we then get that AV is quadratic in A{z). But that means that it 
must be the same as its expansion to quadratic order in A, and so must be the same as the 
result in the linearized theory, giving 

AV[A{z)]^ f \A\^,,Al (A6) 

In k space, one may therefore write 

AV[A{z)] = I \ A^ky n,,(0, k) A'^ik). (A7) 

Since the Fourier transform oi A — A{z) has support only for fc's proportional to e^, we can 
replace 11*^ (0, k) by the matrix of constants 11^(0, e^). The effective potential for A — A{z) 
is then local in x and may be written in the form (3.3). 

2. For SU(2) with k± = 0, ^ and A^ = 

In this section, we will simplify the potential V of (3.5) under the assumptions that (i) A 
lies in an SU(2) subgroup of color SU(3) and (ii) A^ = 0. First, let's consider the restriction 
to an SU(2) color subgroup. Then /"^^ can be replaced by e"^'^, giving 

Vsu(2) = l/[(A'A')'-A^A^4]-i/^'(^S^S + ^^^^) 

= i/|[tr(>lM)]2-tr[(>lM)']}-i/x'tr(A^p(^^)A), (A8) 
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where A is the 3x3 matrix of Aj, 



[Al Al Al] 

\Al Al AlJ 



(A9) 



p{xy) jg ^j^g projection operator onto the xy subspace of directions i, 




(AlO) 



The potential is symmetric under (i) spatial rotations in the xy plane, and (ii) color rotations. 
These symmetries can be summarized as 



where C is a color rotation for the adjoint representation of SU(2), represented by any 
3x3 real, orthogonal matrix with dct C = 1, and i? is a 3 x 3 real, orthogonal matrix 
representing a rotation in the xy plane. By a color transformation, we can assume without 
loss of generality that the matrix A is symmetric.^^ Now make the remaining simplifying 
restriction that A^ — 0. Then the symmetric Aia is zero except for i — x,y and a — 1,2. We 
can then diagonalizc the symmetric A in this 2x2 subspace by a simultaneous space/color 
rotation of the form A RAR^ . So, without further loss of generality, we may write 



This gives (3.7). 

APPENDIX B: UNSTABLE, NON-ABELIAN, STATIC WAVES 

In this appendix, we show the existence of non-linear, non-Abelian, static wave solutions 
to the effective theory (Al) of the soft modes. These solutions are classically unstable, 
however, to the Abelianization discussed in this paper. In the zero momentum limit (/c — > 0), 
we will see that they correspond to the saddle points in Fig. 2. 

The solutions are simple generalizations of propagating non-linear wave solutions found 
for equihbrium plasmas by Blaizot and lancu [29]. We will look for solutions of the form 
A — A{K • x), for some 4- vector K. [We will shortly speciahze to K in the z direction, so 
that A = A{z).] As discussed in Appendix A, Blaizot and lancu then showed that, in this 
case, the only effect of hard particles at leading order is to induce the HTL self-energy for 
the soft fields. The soft equation of motion is then 



Using singular value decomposition, we can write A = DR for any A, where L and R are orthogonal 
matrices and £> is a diagonal matrix. Then AR^ L = DL is symmetric. So, a color rotation of A by 



A RAC^, 



(All) 




(A12) 



(Bl) 



R^L makes A symmetric. 
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Following one of the possibilities studied by Blaizot and lancu, let's focus on an SU(2) 
subgroup of the gauge group and further restrict attention to gauge fields of the form 

^^(n • x) - [e^ai + e1^-^6a2]h(K ■ x), (B2) 

where e|^^ 2) ^ basis of spatial polarizations orthogonal to k. [This is essentially just (A12) 
with 4>i = 4>2 = h{z), if one chooses K along the z axis.] The equation of motion (Bl) then 
becomes 

{u^ - k^)h" + Ut{ou, p)h + g^h^ = 0, (B3) 

where h" indicates the second derivative of h with respect to its argument K ■ x. 

Blaizot and lancu considered cq. (B3) for positive, which is the situation for prop- 
agating modes in equilibrium situations. There are no non-trivial static solutions {uj = 0) 
in this case because 11^(0, fe) = in equilibrium. The difference in this paper is that we 
are interested in cases where 11^(0, k) can be negative because of the Weibel instabihty for 
anisotropic distributions. As a result, static solutions can exist to (B3). Taking k to be in 
the z direction, and writing 11^(0, k) — —jj? < 0, the static {uj — 0) case of (B3) is 

-p,-^'h + g'h'^Q. (B4) 

The solutions are 

for any m with < m < 1. Here sn.{u\m) is the Jacobi elliptic sine of argument u and 
parameter m. 

Now note that (B4) is just the equation for static solutions of a scalar theory with potential 
V{h) = —\ii'^h? + \g'^h'^. This double-well potential is simply the potential of Fig. 2 along 
the diagonal 01 = 02- The minima of V{h) correspond to the saddle points of 1^(01, 02)- In 
the ansatz (B2), corresponding to (f)i — (f)2, we simply cannot see the unstable directions of 
the saddle point. 

For m <C 1 in (B5), the solution has small- amplitude variations around h = as 
one changes z. As m — > 1, the solution becomes a regular array of far-separated kinks 
and anti-kinks, where the standard kink solution, centered on the origin, is h{z) — 
(n/g) tanh(/x2;/\/2). In this case, the kinks interpolate back and forth diagonally between 
the saddle points of Fig. 2, and in most of space, the field is nearly constant near one saddle 
point or the other. And so the constant-field instabilities discussed in the main text will lo- 
cally afflict these solutions. For smaller m, where the kinks and anti-kinks begin to overlap, 
the potential energy can be further minimized by annihilating pairs of kinks and anti-kinks 
(that is, by dimerization) : so again, these represent unstable solutions. 



APPENDIX C: CAN SOFT EXCITATIONS STOP INSTABILITY GROWTH? 

In this appendix, we give some qualitative arguments against the concerns raised in 
section IIIB. We will argue that one should not expect effective masses generated by the 
soft sector to be able to prevent the growth of the system down along the valleys of the 
effective potential shown in Fig. 2. Because we discuss theories in different dimensions in 
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FIG. 15: Like Fig. 2 but including a hypothetical contribution from soft modes that cancels the 
due to hard particles. This potential is zero along the axis. Qualitatively, the creation of 
just slightly more power in soft modes would cause it to curve up instead. 

this paper, we shall keep the spatial dimension d general. For simplicity, we will restrict our 
discussion to the purely scalar theory of our warm-up toy model (3.9). 

Imagine that the system starts with very tiny random magnetic noise, and the modes 
with k < fi then grow due to the instability. In a time of order fi~^ (neglecting factors of 
logarithms), the instabilities grow to the size of the soft non-perturbative scale fi/g. Let's 
start by supposing that the growth of the system did hang up at this scale. That is, suppose 
that the excitation of soft modes has built up to a point such that their contribution to 
the effective potential now cancels the destabilizing — quadratic coefficient contributed 
by the hard modes, so that, crudely speaking, the total effective potential at this time looks 
somewhat like Fig. 15 rather than Fig. 2 (which accounted only for the effects of the hard 
particles). The potential energy density released at this point would be of order fi^/g'^, 
which will also be the kinetic energy density of the soft field. At least initially, this energy 
would be in Fourier components with ~ /i and with non-perturbatively large occupation 
numbers of order l/g"^. If the system were hung up, however, then it would quickly start to 
equilibrate by the strong non-perturbative self-interactions of these modes. 

If the hung-up system were equilibrated, then the kinetic energy density fi'^/ g"^ would be 
j'd+i^ giving T ~ (yU^/(7^)^/*^'^+^\ This gives T ^ fi for weakly-coupled theories jj.'^'^g'^ <^ 1 
in d spatial dimensions. The energy would then be dominated by modes with momentum 
of order T ^ fx. The drive toward thermalization would therefore push the energy of the 
system from modes with /c ~ /i to modes with higher and higher momenta, until the system 
eventually thermalized (if it in fact remained hung up at ~ yu/gf). 

Let A be the momentum scale that dominates the total kinetic energy at any time, which 
starts at order /i. Due to the drive of the system toward thermalization, A would grow 
several times in a time of order Now consider a time when A 3> /x. At this time, 

the kinetic energy density [i^ j g^ ^ which must be of order is distributed among modes 
with momentum ~ A. So now ~ [i} j gk. By mean field theory, the corresponding 
contribution to the squared effective mass is 

(Am2)A ~ g^^ ~ ^. (CI) 

Were A ~ /i, this might or might not be enough to compensate for the unstable in 
the potential. However, since thermalization drives A larger and larger, one will soon have 
-|- (Am^)A < 0, and so instability to further growth, in a time scale of order . 
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There was an implicit assumption above that the momentum scale A which dominates the 
kinetic energy is also the scale giving the dominant contribution to Am^. In low dimensions, 
this need not be the case. But even if all the energy were in modes k ^ /i, corresponding to 
A ~ // above, we have seen that the resulting (Am^)^ is order //^. Once most of the energy 
has left modes // ~ /c, then (Am^)^ will be <^ /i^. So neither (Am^)^ nor (Am^)A will be 
enough to stop the continued growth of the instability.^'^ 



APPENDIX D: LATTICE VERSION OF TOY MODEL 



The lattice version of our 1+1 dimensional toy model, with spatial lattice spacing a, 
consists of the following elements. There are adjoint scalar fields (p^ and (py, corresponding 
to Ax and Ay, which live on the sites of the spatial lattice, represented as NxN traceless 
Hermitian matrices of su(A^). There arc SU(A^) group elements U living on the spatial 
links, corresponding to exp{iaAz) and represented as NxN unitary matrices. Then there 
are the conjugate momenta Ux and 11^ (corresponding to A^ = —E^ and Ay = —Ey) and 
E (corresponding to E^), which are all represented as NxN traceless Hermitian matrices 
living on the temporal links. In this appendix, we take g — 1 and ji — 1. For time step e, 
our evolution equations are 

(l>a,s{t + e) = (t>c,,s{t) + en„,,(i + |e), (Dl) 
U Mt + ^) = U Ut) expf-ieaEsit + le)), (D2) 



\ V s+2 ^"'"2 *~2 «-2^ 

+IJ''^(l>a,s - [(f>a,s, [(pa,s, <Pa,s]]j , (D3) 



where s is an integer site index, a runs over x and y, and a denotes y and x, respectively. 
Gauss's Law is 

a-^(Es{t - \e) - U^_i{t) Es-i{t - |e) U^^_i{t)) +J2'[(PaAt) , ^aAt ' b)] = 0- (D5) 



Here's a different form of the analysis. Consider a system that has equihbrated for modes fc < A <C T, so 
that the occupation numbers are f{k) ~ T/Ek for fc < A and are small for larger k. Then the kinetic energy 
density is order J d'^k Ek,f{k) ^ TA'^, which we set to /i^/g^ as above. From this, we obtain T ~ ji'^/g^A'^. 
The contribution to the effective mass from ah the soft modes is Am^ ^ 9^ J d'^kf{k)/Ek ~ g^T J cfik/El. 
Using the previous expression for T, this is Am? {fi'^/A'^) J d'^k/E'^. For d > 2, the UV contribution 
dominates, giving Am^ ^ fi^/A"^. For d < 2, the IR contribution dominates, giving Am^ ^ fi 
d = 2, Am^ (m^/A^) log(A//i). In aU of these cases, Am^ <C iJ? when A > /i. 
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The approximately conserved energy (which is exactly conserved for e — > 0) is 
£ ^aY^Wsl) + ^tr(n^,J + ^ (('^".^ - 

s '- a a ^ 2 / 

- P^H<t>l,s + <t>ls) + tr((i[0.,„ <t>y,s\y) 



■ (D6) 



In our simulations, the total energy of the system is very close to zero, because we start 
with very small initial fields. One can get a measure of how well energy is conserved in 
our simulations by comparing the change £{t) — £^(0) in the energy with time to the size 
of the potential energy V{t) given by the last three terms of (D6). By this measure, the 
simulations presented in this paper conserve energy at the level of 0.1%. We have checked 
that reducing the time step e from the value e = 0.001 used in our simulations (i) does not 
significantly change any of our simulation results, and (ii) indeed further improves energy 
conservation. 

Our initial conditions (4.7) are 

i(0) = l, 0a,.(O) = O, i?,(-ie) = 0, (D7) 

s-l-2 

{K,s{-¥) ^^A-¥)) = ^ ^"'^apS^s', (D8) 

(Jj 

where {C,^) = A^/a means choose ^ to be a Gaussian random number centered on zero with 
standard deviation a — A/a^/^. In our simulations, a — 0.01. The initial energy density po 
is of order 

Po ~ 0^ ~ cr^ (D9) 

and is dominated by stable UV modes. 

We should check that our simulations will not run into a common difficulty in naively 
applying classical evolution equations to many problems, which is exemplified by the ultravi- 
olet catastrophe of the classical blackbody spectrum. In thermal equilibrium, classical field 
theory (unlike its quantum counterpart), populates arbitrarily high momentum modes with 
real excitations. The interaction of these excitations with low-momentum modes of interest 
can cause arbitrarily large effects, in the continuum limit, on effective masses, damping rates, 
and other properties of the soft fields. (See, for example, Refs. [52, 53].) Our simple white 
noise spectrum (D8) populates arbitrarily high-momentum modes, so one might worry that 
there could be a similar problem here.^^ As an example, let's check the effective masses 
of soft modes due to interacting with UV modes of our lattice. Such masses arise from 
diagrams like Fig. 3. There arc a number of ways to compute their scale. Treating the UV 
lattice modes with kinetic theory, for example, 

Am^^g^ f ^. (DIO) 
Jp P 

where the integration is restricted to the UV modes whose effects are of interest. Previously 
in this paper, fp has referred to hard particles. In this context, it refers instead to the 
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If there were, it could be fixed, of course, by cutting off the spectrum of the initial fluctuations. 
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UV modes of our classical continuum effective field theory of the "soft" modes. The energy 
density of these UV modes is p ~ Epfp. Parametrically, then, 

Am^ ^ ^ ~ (Dll) 

The effects on the soft physics of interest will be small if Am^ <^ /x^. Adopting our conven- 
tions g — 1 and a* = 1, and using (D9), this condition becomes 

a^a^ < 1. (D12) 

In our simulations, cr^a^ = 2.5 x 10~^. 



APPENDIX E: ULTRA-RELATIVISTIC BENNETT SELF-PINCH 

Traditional plasma physics is often complicated by the presence of multiple scales arising 
from the fact that ions are much heavier than electrons. Such complications vanish for 
electron-positron plasmas, or quark-gluon plasmas. In this appendix, we will summarize 
the self-pinching solutions to the coUisionless Vlasov equations found by Bennett [47] in the 
simplifying situation where opposite charges have the same mass. This includes the ultra- 
relativistic limit, where particle masses are ignored. For simplicity of notation, we will focus 
on Abelian gauge theory. A class of self-pinching solutions for a plasma of charges ±e and 
mass m is given by^^ 

f±{x,p)^p{r^)F±{p) (Ela) 
in cylindrical coordinates {r±, 0, z), where F±{p) are boosted thermal distributions 

(with Zm Zq — T^/tt^ in the ultra- relativistic limit T ^ m), and p{r±) is the radial 
profile 

The constants po, and u are arbitrary. One may check explicitly that (El) satisfies the 
time-independent coUisionless Vlasov equations 

V ■ V/± ±e{E + vxB)- Vpf± = 0, (E2) 

V-E^e /(/+-/-), (E3) 
Jp 

VxB^e f v{U-f.), (E4) 



p 



One may derive this equation by starting from the ansatz (Ela) for some unknown funetion p{r±), then 
use Ampere's Law to determine B, and then use the coUisionless Boltzmann equation for f± to obtain a 



differential equation for p. Solving that equation yields (Elc). See Ref. [47]. 
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with E — and 

4Tbpor_L 
e7„n(l + 6po^l) 

This solution represents a finite-width beam of positively charged particles moving in the 
+z direction and an exactly similar beam of negatively charged particles moving in the —z 
direction. At their centers, each of these beams has density po and looks like a thermal 
distribution with temperature T boosted to a net beam velocity of magnitude u. The beam 
velocity u can take any value < m < 1, even in the ultra-relativistic limit, since u is the 
average velocity of particles in the beam. 
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